Static and dynamic lengthscales in a simple glassy plaquette model 
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We study static and dynamic spatial correlations in a two-dimensional spin model with four- 
body plaquette interactions and standard Glauber dynamics by means of analytic arguments and 
Monte Carlo simulations. We study in detail the dynamical behaviour which becomes glassy at 
low temperatures, due to the emergence of effective kinetic constraints in a dual representation 
where spins are mapped to plaquette variables. We study the interplay between non-trivial static 
correlations of the spins and the dynamic 'four-point' correlations usually studied in the context of 
supercooled liquids. We show that slow dynamics is spatially heterogeneous due to the presence of 
diverging lengthscales and scaling, as is also found in kinetically constrained models. This analogy 
is illustrated by a comparative study of a froth model where the kinetic constraints are imposed. 
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I. INTRODUCTION 

Recently there has been considerable interest in the 
extent to which the slow dynamics of glass-forming liq- 
uids may be understood as a result of constraints on 
the kinetic properties of a system, rather than on its 
thermodynamics 0, 0, 0. It has been shown that ki- 
netically constrained models with trivial thermodynamic 
properties |J, |3, |6( show a slowing down at low temper- 
ature, accompanied by the stretched exponential relax- 
ation and dynamically heterogeneous behaviour charac- 
teristic of glass-formers. The diverging timescales arise 
from dynamical fixed points at zero temperature 0, , 
associated with diverging dynamical lengthscales. 

However the only degrees of freedom in these models 
are phenomenological 'mobility fields' Excitations in 
the mobility field represent regions of the glass-former 
where motion is possible. Presumably the energy bar- 
riers that prevent relaxation in the system are smaller 
than average in these regions, but there remain questions 
as to how these degrees of freedom are related to other 
physical properties of the glass-forming systems. In par- 
ticular, is local mobility associated with some (possibly 
complicated) local static ordering? Is it always neces- 
sary to consider dynamic correlators? Is a mapping from 
interacting particles to a mobility field feasible? 

An interesting class of models in which mobility fields 
emerge naturally from more familiar interacting degrees 
of freedom was investigated in Refs. 0, |n| . In this par- 
ticular case, mobile regions are associated with defects 
in the spin fields and these defects have a density that 
vanishes at zero temperature, accompanied by a diver- 
gence in the relaxation time of the system. These models 
therefore provide a specific, yet informative example of 
effective kinetic constraints arising from simple dynami- 
cal rules in conjunction with familiar Hamiltonians with 
multi-spin interactions with no quenched disorder. 

In this paper we investigate correlations in the two- 
dimensional plaquette model 0, 0, 0, ^| . The ther- 



modynamic properties of the model are those of non- 
interacting point-like excitations. However, when these 
energetic properties are combined with a simple spin-flip 
dynamics, the dynamical evolution at low temperatures 
is controlled by large energy barriers. The relaxation 
times diverge in an Arrhenius manner consistent with 
the behaviour of strong glasses 01 ■ This is accompanied 
by a diverging lengthscale associated with dynamic cor- 
relations. Our study therefore illustrates how a mobility 
field might arise from a physical interacting spin system. 
The resulting dynamic behaviour is very similar to the 
physics of kinetically constrained models as we show by 
comparatively studying a two-dimensional kinetic model 
inspired by the physics of covalent froths |l5j . The major 
difference is that in the latter case dynamics is directly 
expressed in terms of the fundamental excitations and 
kinetic constraints are therefore imposed by hand. We 
show that the dynamical correlations in both models can 
be explained in terms of freely diffusing excitations. 

In section[n]we define the models, study static correla- 
tions and the representation in terms of diffusing excita- 
tions. In Section IlIII we study dynamic correlations and 
Section ITVl contains our conclusions. 



II. PLAQUETTE AND FROTH MODELS 
A. The plaquette model 

The Hamiltonian of the plaquette model reads 0, 111) 
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(1) 



where the {o~ij} are Ising spins on a two-dimensional 
square lattice with free boundary conditions. We spec- 
ify simple spin-flip dynamics with rates given by Glauber 
probabilities. We denote the linear size of the lattice by 
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L, so that there are N = L 2 spins in the system. To 
demonstrate that the thermodynamic properties of this 
model are trivial, we make the change to plaquette vari- 
ables, {Pij}, defined as 



Pij — &ij&i+l,j&i,j + l<?i+l,j + l- 



(2) 



The plaquette variables form a dual representation of the 
spin system. For a system with N spins, the partition 
sum over these N variables is equivalent to summing over 
(L — l) 2 independent plaquettes, with a final summation 
over (2L — 1) spins that lie on two orthogonal boundaries 
of the system. 

The result is that excited plaquettes, that is those with 
Pij = — 1, are uncorrelated in space because the Hamil- 
tonian becomes trivial in the plaquette representation, 
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H = - 



2 ^ 



(3) 



Excitations have a density [eP + where /? is the in- 
verse temperature. Also, the free summation over bound- 
ary spins means that all states have a degeneracy of 
2 1 . This arises from the symmetry of the Hamilto- 
nian under flipping all of the spins in any row or column 
of the square lattice. A result of the symmetry is that 
only correlation functions that are invariant under this 
symmetry can take finite values, since there is no spon- 
taneous symmetry breaking in this model. For example, 
the two-point correlator 
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ij<~> i+x.j+yl 



9,0 



(4) 



It is clear that this correlator must vanish if either x or 
y is non-zero, since flipping either row i + x or column 
j + y changes the sign of the correlation function without 
changing the energy of the state. 

Interestingly the plaquette model also describes Ha 
the paramagnetic phase of the eight vertex model |17j . 
The additional presence of ordered states in the eight ver- 
tex model makes it a suitable analogy for a glass-former 
with a thermodynamic melting transition. It turns out 
that the plaquette model remains a good description for 
the 'supercooled' states below T c , see Ref. M 



B. Static lengthscales 

We now discuss the aspects of the static properties of 
the plaquette model that are relevant to the physics of 
supercooled liquids and kinetically constrained models. 
We mentioned above that the concentration of excited 
plaquettes in the model is (eP + l) -1 . The interesting 
behaviour in the model occurs at T < 1 (so j3 > 1), in 
which these excitations arc dilute. We therefore define 



« 1. 



(5) 



Two typical configurations of the spins at two different 
low temperatures are shown in Fig. ^ It is clear that 




FIG. 1: Typical spin configurations at /3 = 5 (left) and j3 = 7 
(right). At the lower temperature we identify the excited 
plaquettes with circles, which appear as 'corners' in the spin 
field. At the higher temperature we only highlight some of the 
defects because they are more numerous. The smallest visible 
lengthscale is the lattice spacing. All two-point correlations 
vanish, even between adjacent sites. 



even though all two-point functions vanish by the sym- 
metry discussed above, the system is not in a typical 
paramagnetic state. The axes of the underlying lattice 
are apparent, and their influence is felt even at relatively 
large lengthscales. 

To understand static correlations we note that excited 
plaquettes are topological defects, in the sense that while 
their energy cost is localised, removing a single excited 
plaquette requires changes in the spin field over large 
distances. This is apparent from the behaviour of the 
four-point correlation function 



^fstat 



(x,y) 



ij u i+x,jaij+yVi+ x ,j+y/i 



a, 



(6) 



which is independent of i and j in a translationally invari- 
ant system (or sufficiently far away from the free bound- 
aries of a finite one). We have added the superscript 'stat' 
to emphasize that we consider here static correlations, as 
opposed to the 'four-point' dynamic correlators studied 
in Section UTTl The symmetries of the Hamiltonian mean 
that © is the only non-trivial four-point function that 
can have a finite value. Symmetries also ensure the van- 
ishing of all the disconnected parts of that function. 

Now, the value of the combination 
VijVi+x,jVi,j+yVi+x,j+y is given by the parity of 
the number of excited plaquettes in the rectangular 
region defined by the four spins. To see this, we may 
take x, y > without loss of generality, and write 



i j U i + x , j u i , j + y u i + x , j + y 



x — 1 y—X 

n n Pi+*',i+v' ( 7 ) 

x'=0 y'=0 



where the p^ were defined in (J2J. This rewriting relies 
only on the fact that the spins are Ising variables so that 
afj = 1. The product is over all plaquettes in the rect- 
angular region and is +1 if there is an even number of 
excited plaquettes and —1 if there is an odd number. 
Thus, if there is only one excited plaquette in the rect- 
angle, then the product of the four spins is equal to — 1, 
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regardless of x and y. This is a signature of a topological 
defect. 

When excited plaquettes are dilute the number of exci- 
tations in a rectangular region of area \xy\ is a Poissonian 
random variable with mean c|a;?/|. The probability that 
this variable takes an even value is (1 + e~ 2c ™)/2. The 
four-point static correlator is therefore 

Cf at (x,y) = e - 2c ^l, (8) 

from which one can identify the static correlation length- 
scale £| tat = cr 1 / 2 = e^l" 1 . Physically this lengthscale 
measures the typical spacing between excited plaquettes. 
Clearly this lengthscale diverges at T = where the con- 
centration of defects vanishes. 

It is interesting to note that in a standard kinetically 
constrained model where the Hamiltonian is directly ex- 
pressed in terms of the relevant excitations, the distance 
between defects is trivially related to their concentration. 
In the present case, the study of simple structural cor- 
relators, e.g. two-point structure factors, does not pro- 
vide information on these relevant lengthscales; they are 
only revealed when higher-order correlators are consid- 
ered. More physically this means that the typical size 
of black and white domains in Fig. ^ does not depend 
on temperature while the concentration of circles does, 
so that by naively looking at the interacting spins one 
would get the wrong impression that there is no diverg- 
ing lengthscale in the system. Taking the analogy with 
liquids seriously this suggests that while the structure 
factor of supercooled liquids shows no particular trend 
when temperature is decreased, higher-order correlators 
related to some yet unknown local ordering could reveal 
the presence of increasing lengthscales, as assumed in 
several scenarii of the glass transition [||| • 

The problem with the correlator 10 is that it is highly 
specific to the model under study, which has special sym- 
metries. It would be more useful to define static quanti- 
ties which are independent of the model under study but 
still reveal the presence of growing lengthscales. An idea 
is that in the presence of diverging lengthscales, large 
fluctuations can also be expected. Let us consider the 
two-point quantity, 

1 L 

d$ &t (x,y) = j^J2 ^3^+*,]+y ( 9 ) 

where the sum is over a finite region of an infinite spin 
system. The number of terms in the sum is N = L 2 ; 
we will later take the limit of large L so as to extract a 
well-defined measure of the size of fluctuations. 

The expectation value of c| tat is a sum of two point 
correlation functions. Eq. (0J therefore constrains it to 
vanish in the thermodynamic limit (except in the trivial 
case, x = y = 0). In other words, (c| tat (x, y)) = S Xi oS y fi- 
However, even when the expectation value of c| tat is zero, 
its fluctuations are not. One can therefore define the 
following susceptibility, 

xT'i^y) = N[{cr\x,y) 2 ) - (cr\x,y)) 2 }, (10) 



where the disconnected part in fact vanishes by symme- 
try [except at x — y — where it exactly cancels the 
connected part, so that xl tat (0,0) = 0]. The connected 
part is a sum of four point static expectation values, but 
symmetries of the model constrain most of these to be 
zero as well. The only non-zero terms are either of the 
form of or else trivial (for example (cr 2 jCr 2 j,)). As- 
suming that at least one of x and y is non-zero, we write 
out the sums and collect the non-vanishing terms, arriv- 
ing at 

I L 

X S 2 at ( X >V) — Jr X] (VijVi+xj+yCTi'j'Vi+xJ'+y) 
i,i'j',i'=l 

= (1 - S Xt0 )(l - Syfi) + 

1 L 

£ [<koCf at (*-m,y)+ 

<5 y , C<f at (x,j-m)] (11) 

We note that the first term gives a trivial value of unity 
for xl tat ( x J?/) unless at least one of x or y is equal to 
zero. This arises from the four point correlators with 
i = i' and j = f; we will see shortly that x| tat (x, 0) also 
approaches this trivial value in when x becomes much 
longer than the relevant correlation length. 

The definition of c| tat contains a dependence on the 
cluster size, L. However, the physically relevant limit is 
when L is much greater than all correlation lengths in the 
system; in this limit then x| tat (x, u) converges to a finite 
value. From JSJ), we can identify the sums as geometric 
series: taking care not to double count the terms with 
i = i' or j = j' , we find that 

oo 

lim x s 2 tat (x,0) = V e - 2c l^l =coth(c|x|), (12) 

y=-oc 

for x ^ 0. From these fluctuations we identify a second 
static lengthscale 

^ tat = c- 1 = eP, (13) 

which also diverges at T = 0. Physically this length- 
scale represents the mean distance between adjacent ex- 
cited plaquettes in the same row or column of the square 
lattice, and is therefore much greater than the typical 
spacing between excited plaquettes, £| tat . 

In deriving Eq. (|13l) we again made use of the special 
symmetries of the model, which is an undesired feature 
of this treatment. This can be cured by noting that £| tat 
may be accessed through a circular average, 

X s 2 tat (r) = / — X stat (z, y) = — coth(cr) + 1 (14) 

where the integral over 8 is to be understood as a sum 
over the integrand at values of x, y such that r — i < 
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FIG. 2: Simulation data confirming the analytic results. We 
show xl tat f° r V = at three different temperatures (filled 
symbols); the solid lines are the result of 11 2t . We also show 
that xl tat takes the trivial value of unity when neither x nor 
y is zero. We use a cluster of 10 6 spins to avoid finite size 
effects. 



\Jx 2 + y 2 < r + |. Thus the lengthscale £| tat can in- 
deed be extracted from fluctuations of two-point mea- 
surements that do not depend on knowing the exact ori- 
entation of the lattice axes or the special symmetries of 
the model. 

To verify these exact results, we performed static 
Monte Carlo simulations to obtain xf at { x iV) both for 
y = and for the trivial case i,y / where the suscep- 
tibility takes the value unity. The results are shown in 
figure El 



C. Inherent structures and 'droplets' 

We now discuss the above findings in relation to al- 
ternative approaches to the glass transition that involve 
growing lengthscales of a static nature. A problem with 
these approaches is that there exist no numerical or ex- 
perimental indications for the presence of such static 
correlations. However two recent papers have discussed 
methods to detect growing static lengthscales in super- 
cooled liquids that we now discuss in the context of the 
plaquette model. 

Bertin [2l| suggested to consider the spatial structure 
of inherent structures p3 | and studied a one-dimensional 
disordered Potts model to illustrate his ideas. The pro- 
cedure is as follows. Take an equilibrium configuration of 
the system at finite temperature T and quench it to its 
inherent structure. Now fix the orientation of the spins 
on the boundary of a finite portion of the system of linear 
size £. Finally minimize the energy of the finite size sys- 
tem given these boundary conditions. Bertin finds that 
there exists a well-defined lengthscale, £*(T), which is 
such that it is typically possible to find a lower energy 
structure when £ > £*, while the system is already at its 
ground state for smaller sizes, £ < £* . 



In the plaquette model, the effect of quenching to 
T = is that dimers freely diffuse (see below) until they 
get absorbed, and then all motion stops. The inherent 
structure is therefore a structure with single defects only. 
On quenching from equilibrium these defects will be ran- 
domly distributed in space, except that they are never on 
adjacent sites (since that would be a dimer). So, consider 
a rectangular region of the sample of linear size £. We 
may ask, what is the probability that this region min- 
imises the energy, subject to its boundaries remaining 
fixed? Fixing the 4(£ — 1) spins that lie on the bound- 
aries constrains the parity of the number of defects in 
each row and in each column. Thus the (possibly degen- 
erate) ground state of the square region contains at most 
one defect in each row or column, according to their par- 
ities. It is then obvious that when £ > £* ~ £| tat there 
will typically be more than one defect in each line and 
column and the inherent structure will not coincide with 
the ground state of the finite size system. We conclude 
therefore that Bertin's method would successfully deter- 
mine the static correlation lengthscale £if at discussed in 
the previous section. The main difference between our 
model and that of Bertin is that we cannot identify a 
tiling of inherent structures with regions that minimize 
the energy because our defects are point-like objects that 
obviously can not delimit a particular area. 

In a real space description of the random first order 
transition it is imagined that above its static Kauzmann 
transition the system is composed of an assembly of 'en- 
tropic droplets' [Isj. Revisiting this idea Bouchaud and 
Biroli proposed to associate a temperature dependent 
lengthscale, £* (T) , to these droplets in the following man- 
ner p3|. Consider an equilibrium configuration of the 
system at temperature T. Then consider a finite portion 
of the system of linear size £ and freeze its boundary con- 
ditions. Now let this finite size system evolve at temper- 
ature T with its boundary conditions fixed. If £ < £* the 
system should effectively be non-ergodic while dynamic 
correlation functions would go to zero at large times for 
£ > C- 

In the plaquette model this procedure was studied from 
a different perspective in Ref. where the thermody- 
namics of finite size systems was computed exactly for a 
particular choice of boundary conditions. It was found 
that in a finite size system there indeed exists a 'glass 
transition' temperature, which is manifested by a down- 
ward jump in the specific heat which shifts to lower tem- 
perature for larger £. Although this calculation is not 
exactly the procedure described by Bouchaud and Biroli 
where any sort of boundary conditions should be stud- 
ied, it is very close in spirit. Physically, the non-ergodic 
behaviour is due to the fact that when £ is too small the 
system cannot rearrange without altering its boundaries. 
This happens as above when lines or columns only con- 
tain one or no defect, i.e. when £ < ^2 tat and we conclude 
that the procedure of Ref. [2^ would once more yield the 
correct correlation lengthscale, £* ~ £,2 tat - 
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FIG. 3: Sketch of a dimer explaining the origin of the strongly 
anisotropic dynamical correlations in the plaquette model. 
Excited plaquettes are marked with a x sign. Spins are 
marked by small circles but the state of the spins is not shown. 
The dimer diffuses along the x direction with no energy bar- 
rier. Open circles mark spins that can be nipped as part of 
this free diffusive process. Flipping of spins marked by closed 
circles costs at least AE — 2. 



D. Diffusing excitations in the plaquette model 

We now consider the dynamics of the plaquette model. 
The fundamental moves are spin- flips. When a single 
spin is flipped the states of the four plaquettes surround- 
ing that spin must all change. Thus there are five types 
of move, depending on the environment of the relevant 
spin. If a spin has no adjacent excited plaquettes, then 
flipping that spin incurs an energy cost AE = 4, so these 
moves are suppressed by a factor c 4 . Now, if a spin is ad- 
jacent to exactly one excited plaquette, then its flips cost 
AE = 2 and are suppressed only by a factor c 2 . The den- 
sity of such sites is approximately c at low temperatures, 
so these flips are already more significant than those in- 
volving AE — 4. Finally spins adjacent to pairs of ex- 
cited plaquettes can flip without energy cost, AE = 0, 
and that those adjacent to three or four excited plaque- 
ttes have AE < and will therefore relax rapidly to 
the low energy state with one or zero excited plaquettes, 
respectively. 

The key physical point is that excited plaquettes act 
as sources of mobility, since the energetic barriers to spin 
flips are smaller in those regions. This observation al- 
lows us to identify the excited plaquettes as the excita- 
tions in the 'mobility field' by analogy with kinetically 
constrained models. The square lattice of the plaquette 
model has a surprising effect on the diffusion of these mo- 
bility excitations, because pairs of adjacent excited pla- 
quettes diffuse freely, but are confined to one dimension, 
see Fig. |21 As shown in Refs. 0, 0, this is responsi- 
ble for the form of the low temperature divergence of the 
relaxation time, r ~ c~ 3 . 



E. The froth model 

We now introduce a simple model that shares many 
features in its dynamical behaviour with the plaquette 
model. Consider a mobility field that is defined on the 
plaquettes of a hexagonal lattice. Each defect carries an 



FIG. 4: Two possible moves in the froth model. The binary 
variables rn are defined on the hexagons. The + and x signs 
mark sets of four m. A single move involves flipping all four 
of any such set. The 'Feynman diagram' shapes are shown in 
bold. 

energy of unity, so the Hamiltonian takes the trivial form 

1 N 

gfroth = (15) 
i=l 

where there is a defect on site i if rij = — 1; otherwise 
rii = +1. Dynamical moves involve flipping the values of 
four of the n.; on plaquettes that are related by a 'Feyn- 
man diagram' shape, see Fig. 0] This sort of model was 
studied in the context of covalent froths , so we refer 
to it as the 'froth model' in what follows. The similar- 
ity between Hamiltonians (Q and (|15|l is evident but the 
froth model does not have an equivalent spin representa- 
tion. 

Comparisons between the dynamics of the froth model 
and those of the plaquette model are instructive because 
the only difference between their dynamics is the struc- 
ture of the underlying lattice. However, while the dif- 
fusing pairs of excitations in the plaquette model were 
confined to one dimension as a result of the square lat- 
tice, Fig. |3 the equivalent dimers of the froth model are 
free to diffuse across the whole two-dimensional plane, 
Fig.H 

To first illustrate the strong effect of the choice of lat- 
tice on the dynamics we show typical realisations of the 
persistence function j2j of each system in Fig. [5] The per- 
sistence function on site i is defined by bi(t, t') = 1 if the 
state of that site has remained constant for all times be- 
tween t' and t + t', otherwise bi(t, t') = 0. ft is clear from 
Fig- El that the effect of the square lattice is to introduce 
strong anisotropy into the dynamic spatial correlations. 
The remainder of this paper contains discussions of the 
dynamic correlations shown in Fig. [5J 



III. DYNAMIC CORRELATIONS 

A. Dynamic correlations in the plaquette model 

In the previous section we introduced the persistence 
function to quantify dynamics in the plaquette model. 
This was a convenient way to compare dynamics between 
froth and plaquette models. However the presence of the 
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FIG. 5: Realisations of persistence function in the plaquette 
model (top) and the froth model (bottom). Plaquettes that 
have flipped between time zero and time t are coloured black. 
The time t is such that the fraction of black plaquettes is 
approximately 0.4. Both simulations are at /3 — 5 and show 
regions of size 200 x 200. The strong anisotropy in the top 
panel is purely a result of the underlying square lattice in the 
plaquette model. 



underlying spin field in the plaquette model means that 
its dynamics is more naturally investigated in terms of 
the spin-spin autocorrelation 



ij (t,t')=a ij {t , )a ij {t' + t). 



(16) 



The two-time form of this operator makes it more natural 
than the persistence function which depends on the spin 
at all times between t' and t' + t. 

A typical realisation oiaij(t, t') is shown in Fig. EI The 
autocorrelation function retains slightly more informa- 
tion than the persistence function. If the system makes 
an excursion to an excited state before returning to its 
initial one, then the autocorrelation function records the 
fact that the overall state has not changed, but the per- 




FIG. 6: Typical realisation of the autocorrelation in the pla- 
quette model. Sites with dij = — 1 are coloured black. The 
inverse temperature is (3 = 5 and the timescale is such that 
(aij) — 0.6. This figure is from the same time series as for 
the persistence field in Fig. [K] (top) but the time t is longer 
because most of the contributions to the persistence are from 
dimers that are reabsorbed at their emission point. 



sistence function does not. For this reason, the autocor- 
relation data shows the spatial correlations more clearly 
than the persistence data. 

We now investigate the anisotropic dynamic correla- 
tions in the plaquette model more carefully. Consider the 
two-point, two-time (i.e. 'four-point') correlation func- 
tion 



C 2 ,2(x,y,t) = (aij(t,to)ai+x,j+y(t,to)) 

— {a^ (t, to)) (a,i+ x j+y (t, to)) ■ 



(17) 



In equilibrium this function is independent of (i,j,to) 
since the system is invariant under translations in both 
space and time. It is convenient to discuss also nor- 
malised correlations, 



C2,2{x,y,t) 



C2,2{x 1 y,t) 
1 - (MMo)) 2 



(18) 



so that C 2 , 2 (0,0,t) = 1. 

We are more particularly interested in the scaling 
of 6*2,2 at low temperatures. We find that the scal- 
ing contains contributions from diverging length and 
timescales, consistent with the presence of a dynamical 
critical point at T = 0. For example, the on-site autocor- 
relation function (a,j (t, t')) obeys the scaling relation ^1 
(aij(t,to)) = fi(c 3 t). However the spatial scaling is not 
as simple as this temporal relation. It is clear from Fig.0 
that the presence of the underlying square lattice is rel- 
evant to C 2 ,2 even at low temperatures. 

We have investigated the function 6*2,2 in some detail 
using Monte Carlo simulations with a continuous time 
algorithm j^. Since there is a correlation length pro- 
portional to 1 1 c in the static properties of the system, it 
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FIG. 7: Plots of C2,2(x, y, c 3 t) with y = and with y = x, as a 
function of r = y^x 2 + y 2 . Times are such that c 3 t = 0.07 so 
that (a,ij(t,t')) ~ 0.5. The lengthscale for on-axis correlations 
increases with decreasing temperature, but there is no scaling 
for the off-axis correlator. 

is necessary to use rather large system sizes. We use pe- 
riodic boundary conditions, which require systems large 
enough that typical rows and columns of the lattice con- 
tain at least two excited plaquettes. Otherwise there are 
rather strong finite size effects, as observed in Typ- 
ically we use system sizes of 10 6 spins. The typical error 
bars are of the order of the symbol sizes in the figures, ex- 
cept where noted otherwise. Time is measured in Monte 
Carlo sweeps throughout and all distances are in units of 
the lattice spacing. 

In Fig.[7|we show the strong correlations in 6*2,2 along 
the axes of the model, and their rapid decrease away from 
those axes. In order to show data at different tempera- 
tures, we use times t that scale with c~ 3 so that (a^ (t, t')) 
is the same at each temperature. We see that the length- 
scale associated with C2${x, 0,t) grows rapidly at low 
temperatures while the off-axis correlations measured by 
^2,2(2, x, t) are weak and do not depend strongly on f3. 
In the representation of Fig. [ft] C2,2(x, 0, t) measures the 
typical length of the rod-like objects, which gets large at 
small temperatures. However, the absence of correlations 
away from the axes of the lattice indicates a lack of cor- 
relation between those objects. The microscopic origin 
of this behaviour is in the one-dimensional diffusion of 
pairs of excited plaquettes sketched in Fig. 

We now consider the full time and space dependence 
of C2,2(x,0,t) in a little more detail. It is convenient 
to start the discussion by focusing on the corresponding 
dynamical 'four-point' susceptibility pH l2fj| 

X2,2(t)=^C 2 , 2 (a : ,j/,t). (19) 

xy 

This quantity is proportional both to the typical corre- 
lation area and to the strength of the correlations. From 
Fig.[7|we expect C 2 ,2(x,0,t) ps p(t)f[x/&,2(t)] at large x 
(see also below) so that we may write 
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FIG. 8: Top: Plots of X2,2(t) at various temperatures. Mid- 
dle: Rescaling of time to c 2 t showing saturation at c 2 r sa t ~ 1. 
Lines are a power law, t ' , consistent with diffusive motion, 
and a horizontal plateau showing saturation at c\2,2 ~ 0.6. 
Bottom: Rescaling of time to (?t showing decrease of corre- 
lations due to decreasing p which scales in the same way as 
the spin-spin autocorrelation functions. (The autocorrelation 
functions are displayed as the data that collapse on the dotted 
line.) 



In this expression df is the fractal dimension associated 
with the correlations, £2,2^) the linear size of a correlated 
region, £ 2 2 its area and the quantity p(t) measures the 
strength of the correlations. We show X2,2(i) at several 
temperatures in Fig. [S] (top). There are three distinct 
regimes in this function: a power law increase at small 
times; a broad maximum whose width increases with 
decreasing temperature; and a rapid decrease at larger 
times. 



X2,2(t)»p(tK2,2W d '. 



(20) 



These data are easily interpreted in terms of the exci- 
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tations described in Section [H] At short timescales the 
rod-like objects shown in in Fig. [S] grow in a diffusive 
manner. These objects have dt = 1, and we therefore 
expect X2a{t) ~ (t 1 ' 2 ) d f~ 1 , as observed in Fig- El More- 
over this behaviour is independent of c since a pair of 
excited plaquettes can diffuse at no energy cost. 

At larger times, dimers that have diffused can be ab- 
sorbed on encountering an isolated excited plaquette, 
preventing further growth of the rods. The dynamic 
lengthscale therefore saturates to the mean distance be- 
tween isolated defects along one direction of the square 
lattice. This is precisely the physical content of the static 
lengthscale £| tat ~ c _1 discussed above. Since the pro- 
cess is diffusive, saturation of £2,2 to £| tat takes place at 
a timescale given by 

T S at~C- 2 . (21) 

Saturation at time T sa t is observed in the numerical data, 
as shown in the middle frame of Fig. |SJ 

Finally, after saturation there is little change in the 
correlations until the susceptibility starts to decrease 
for times t ~ c~ 3 . Scaling in this late-time regime is 
governed by the time dependence of the factor p(t) in 
Eq. (|2L)[1 which starts to decrease significantly when the 
rods start to overlap. When two rods cross the site at 
which they intersect has flipped twice. It therefore has 
fly = +1, unlike the rest of the rod which has ay = — 1. 
This effect takes place on the timescale set by the au- 
tocorrelation function, which scales as t ~ c~ 3 3> T sat . 
This behaviour is confirmed in the bottom panel of Fig. El 
where the scaling of (a^) is also shown. 

B. Dynamic Scaling in the plaquette model 

We now discuss the spatial structures in C2,2(#, 0, i), 
and their scaling. At small times, pairs of defects diffuse 
in one dimension. Dynamic correlations only depend on 
x/£,2,2(t) ~ x/t 1 / 2 , and the scaling of x is consistent with 
df = 1. In this time regime we may identify a dynamical 
exponent, 2 = 2. At larger times, T sat <C t, the dynamical 
lengthscale saturates at £2,2 (i) ~ t;^^- This saturation 
of the dynamic lengthscale reflects the fact that the typi- 
cal distance that a dimer travels before being absorbed is 
given by £| tat . We may identify this length as a 'mean free 
path' for dimers diffusing along one dimensional paths. 
Thus the structure of dynamical correlations is explained 
in terms of the constrained one dimensional motion of 
diffusing pairs of defects. The behaviour may be sum- 
marised as 

C2 > 2{XM - \ h b {cH)f 2c {cxl cH»l. ( 22 ) 

We illustrate this behaviour in Fig. [5] The function 
/2a is shown in the top panel, along with our theoreti- 
cal prediction for it [Eq. (|23|l . discussed below]. In the 
middle panel, we show the spatial correlations at vari- 
ous times that are all greater than the saturation time. 
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■ P=4.5, (=10 




U.J 

cx 



FIG. 9: Top: Data collapse of 62,2(2:, 0, i) on reseating of 
space by t 1 ' 2 independently of c. The solid line is the pre- 
diction for this function, assuming free (one-dimensional) dif- 
fusion of pairs of excited plaquettes, valid for times t <C r sat 
[Eq. fZft]. Middle: Scaling of C 2 ,2(x,0,t) for t > r sat . We 
plot temperatures between (3 = 4 and (5 = 5.5 for different 
scaled times c 3 t = 1,2,3,4,5 (top to bottom). Errors are 
smaller than the distances between traces. Bottom: Collapse 
of the data of the middle panel using C2,i(x, 0, t)/[cx2,2{t)] 
showing the scaling function f2c(x) in Eq. 1221 

The spatial dependence of the correlations is not chang- 
ing, while their strength gets weaker and scales with the 
relaxation time c -3 . Finally, we may normalise the cor- 
relations by their spatial integral as a way to extract the 
scaling function ji c {cx) in the bottom panel of Fig. EJ 

Note that these scaling plots contain no free parame- 
ters. There remain small deviations from scaling arising 
from the presence of the lattice at small distances and 
from the fact that the two timescales r sat and r are not 
infinitely separated. However Figgis clear evidence that 
the dynamical correlations are well-described by inde- 
pendent random walkers in one dimension, with a mean 
lifetime of approximately r sat . 

We may evaluate exactly the function f2a(x 2 /t), which 
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gives the value of 62,2(2;, 0, f) for times much smaller 
than the saturation time. The four spin function 
O'ij{t,to)ai^. x .j(t,to) takes the value —1 if either (1) a 
dimer starts between the sites (i, j) and (i + x, j) and dif- 
fuses out of that region, or (2) a dimer starts outside that 
region and diffuses into it, or (3) a dimer diffuses along 
the y direction, flipping one of the spins, or (4) some 
more complicated process occurs, involving two dimers 
or the absorption or emission of a dimer by an isolated 
excited plaquette. 

Assuming that the dimers diffuse freely then the first 
three processes mentioned above are easy to analyse by 
manipulation of the diffusion equation in one dimension. 
We arrive at the following result for 6*2,2(2;, 0, t): 



(23) 

where erfc(x) = (2/y/n) dy e~ v . 

In deriving this result we assumed that only dimers 
present at t = to affect the spins, and that these dimers 
diffuse freely throughout the system. This assumption 
breaks down when the dimers present in the initial state 
start to encounter isolated defects, where they are ab- 
sorbed. This is the physics behind the collapse in the 
middle panel of figure and happens when c 2 t reaches a 
value of order unity. The absorption of the dimers cuts 
off the growth of the dynamical lengthscale, leading to 
shorter range correlations than those predicted, as is vis- 
ible from the top panel of Fig. [§] In particular, we note 
that the correspondence between theory and simulation 
is best for small times (compared to r sat ). 

A final comment on dynamic scaling and critical be- 
haviour. Usually near critical points, the behaviour is 
determined only by long lengthscales, the symmetries of 
the Hamiltonian, and dynamically conserved quantities. 
In the plaquette model, the short lengthscales associ- 
ated with the underlying lattice are relevant even at very 
small temperatures, as can been seen by the strongly 
anisotropic behaviour of 62,2(2:, y, t). However the van- 
ishing of the static correlation function, 6f tat (2;,y), away 
from the axes of the underlying lattice is a result of the 
symmetries of the Hamiltonian. The similar reduction of 
62,2(2:, y, t) has its origin in a conserved quantity, namely 
the number of excited plaquettes in every row and col- 
umn of the square lattice which is conserved modulo 2. 
This is the reason for the confinement of dimers to one 
dimension. The absence of any such conserved quantiy 
in the froth model explains why it is so different from the 
plaquette model. 



C. Suppression of diffusion in the plaquette model 

This subtle difference between froth and plaquette 
models can be made more spectacular when an aging 
situation is considered. It is known jl5j that the aging 
behaviour of the froth model after a quench from high 
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FIG. 10: Time decay of the density of defects, u, after a 
quench from infinite temperature. (Top) Data for /3 = 5 and 
/3 = 6. The initial plateau at short times is a signature of the 
onset of activated dynamics, common to many facilitated spin 
models. As the density of excited plaquettes falls, the system 
enters a state in which the vast majority of spins must over- 
come an energy barrier in order to flip. Thus the dynamics 
slows down, see |(|. The final plateau is equilibration at a 
density (1 + e* 3 ) -1 . (Bottom) Quench data from long times, 
at (3 = 5, 6, 7; these collapse as a function of rescaled time c 2 t, 
until such time as the system equilibrates. The dotted line is 
a fit: Ay/\og{BcH)/(c 2 t) with (A, B) = (0.02, 90). 



temperature is consistent with diffusion and annihilation 
of defects, A+A — ► 0, and has therefore an energy density 
that decays as (\ogt)/t in two dimensions. 

In the plaquette model, the conservation of the parity 
of the number of excitations in each row or column im- 
plies that a single defect cannot simply diffuse. Accord- 
ingly, we show in Fig. QJJ that the aging behaviour of the 
plaquette model is moved out of a simple annihilation- 
diffusion regime, and the energy density decays much 
more slowly than (log t)/t. In Ref. [16(, we showed that 
this decay can also be fitted satisfactorily by a power law 
£-0.45 ( ar i g j n g from an decay equation dtu oc u 2 2 where 
we attributed the anomalous exponent to a fluctuation 
correction). We will argue below that the fluctuation 
correction is in fact logarithmic; we also note that the 
same decay was recently fitted on a much smaller time 
window as i -1 / 3 |13|| . invalidated by our simulations per- 
formed on much larger timescales. 
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To explain this behaviour, we recall that in two- 
dimensional annihilation-diffusion processes the effective 
rate equation is 

+ (24) 

where u = (H/N + 1/2) is the density of excited plaque- 
ttes, and D is the diffusion constant. In the plaquette 
model any diffusive step in fact involves movement of 
at least two excited plaquettes. The simplest such hori- 
zontal moves involve isolated excited plaquettes in either 
the same row, or in adjacent rows, as in Fig. where 
one isolated defect has moved by emitting a dimer. This 
pair must then be absorbed by the other isolated excited 
plaquette, or the system will revert to its original configu- 
ration. If the separation of the two isolated excitations is 
d ^> 1 then the probability of the dimer traveling across 
the gap is approximately d _1 01 . The diffusion con- 
stant is therefore proportional to e _2 ^(d _1 ) where the 
exponential prefactor arises from the activation barrier 
to dimer creation. 

In equilibrium one has (d^ 1 ) ~ l/£! tat ~ c - The 
diffusion constant and inverse relaxation time are both 
proportional to c~ 3 . Out of equilibrium we should also 
consider correlations in the density of excited plaquettes. 
However these correlations take place on a lengthscale 
ss m -1 / 2 , which is much smaller than the likely values of 
d. Therefore we neglect these fluctuations to get D cx u. 
Substituting into Eq. I|24|). we find that 

«~V¥- (25) 

This prediction is in reasonable agreement with Fig. ITU1 
However very many decades of time would be necessary 
to confirm beyond any doubts this particular form of log- 
arithmic corrections. 



D. Dynamical correlations in the froth model 



model, the scaling of C' 22 {x, y, t) is very similar to that 
of C 2t2 (x,y,t), compare Figs. and El Therefore it is 
sensible to compare C 2 2 for the froth model with C 2 ,i 
for the plaquette model. Further, since we will find C 2 2 
to be isotropic for all but the smallest lengthscales in the 
froth model, we write its circular average as C 2 2 (r, t) and 
use this latter correlation function in what follows. 

To generalise the arguments of section llll Bl to the froth 
model, we first recall that for random walkers in two 
dimensions, the typical distance traveled in time t scales 
as i 1 / 2 , but the typical number of sites visited by the 
walker scales as t/logt. Logarithmic corrections arise 
because d — 2 is a marginal dimension for free random 
walks, see Ref. |27j and Refs. therein. 

We expect therefore the following generalisation of the 
three regimes observed in the plaquette model. At small 
times we expect the lengthscale to grow as £/ 2 2 ~ * an d 
the susceptibility to grow with the number of visited sites, 
X2.2 ~ t/logt. Saturation occurs at a lengthscale con- 
trolled by the mean free path of diffusing dimers, which 
scales as c -1 / 2 and represents the separation of free ex- 
citations which absorb the dimers. We note that this 
lengthscale is smaller than the mean free path of the pla- 
quette model, which was the typical separation between 
single defects in the same row or column of the square 
lattice. Since the scaling is diffusive, saturation will occur 
at a time, r s ' at / log r s ' at ~ c _1 . The decay of the suscep- 
tibility will finally be controlled by the decay of the per- 
sistence function on a timescale given by r/logr ~ cr 2 . 

To illustrate these predictions, we first show the per- 
sistence function in Fig. ^2 plotted as a function of the 
naively rescaled variable c 2 i, and also against the more 
appropriate variable c 2 t/\ogt. Since we work over only 
a few decades of time, we may also obtain reasonable 
fit to a form c t, but the logarithmically corrected vari- 
able which has no free parameters appears indeed as the 
appropriate one. 

In Fig. ^] we show the susceptibility associated with 
the persistence function 



Having described the dynamical correlations of the pla- 
quette model in some detail, we can understand the cor- 
relations in the froth model rather easily. There is no 
restriction on defect concentration in rows and columns 
of the hexagonal lattice, so the anisotropic behaviour of 
the plaquette model is absent. However the phenomenol- 
ogy is a simple generalisation of the behaviour described 
above. 

In the absence of spin degrees of freedom, we consider 
the persistence function bij(t,t') instead of the autocor- 
relation function as a local dynamical function. That is, 
we define 



(bij(t,to)bi+ x j+ y (t,to)) — (bij(t,to)y 
(MMo)) - (MM0)) 2 



Ci 2 {x,y,t) 

\ '■' 1 1 \ < - m ! ; / \ " / / v. ' ■ ' 1 ) ; / 

where the prime indicates that we deal with persistence 
rather than autocorrelation functions. In the plaquette 



X2,2(<)=E C 2, 2 (^^)- 

xy 



(27) 



Saturation takes place at a value of x'2 2 that is propor- 
tional to c -1 and the approach to saturation is a univer- 
sal function of ct/ dog t. The susceptibility then decays 
on the timescale of the persistence function. 

The collapse of the persistence and of the susceptibility 
with a single rescaled time variable is consistent with 
their dependence on the number of sites visited by each 
random walker. When considering spatial correlations, 
the situation is slightly more complicated. While the 
relevant areas scale as t/logt, the distances traveled by 
walkers are distributed as universal functions of r 2 /t. We 
therefore expect that while the susceptibility varies as a 
scaling function of t/ logt, the moments of the probability 
distribution should have a normal diffusive behaviour. 
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FIG. 11: Two scaling plots of the persistence function in the 
froth model, showing the effect of logarithmic corrections in 
two dimensions 



That is, we expect 

X2,a (*) = / Zirrdr C' 2a (r, t) ~ tj log t, (28) 



but simultaneously 

/ 2irrdrr n C^ 2 (r, t) 
J 2irrdr C' 2 2 (r, t) 



7/2 



(29) 



These relations were proved recently in Ref. [23 for non- 
interacting random walkers in two dimensions and un- 
normalized dynamic spatial correlators. An alternative 
statement of Eq. I|29|l is that 



C 2 , 2 M) = - h [ — 



t < Tg; 



(30) 



which should be valid for r 1. This scaling relation 
must break down at small distance since both C 2 2 (0, t) 
and x'2 2W are scaling functions of tj logi, incompatible 
with llSfl . 

We test the collapse of Eq. (J201 in Fig. El There are 
tiny deviations at small r 2 /t as expected. Moreover we 
get reasonable agreement with the analytical results of 
Ref. [23 which predict 



2 ^r^(j-i)^' j,2D ' <3i) 
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FIG. 12: Top: Plots of X2,2(0 a t various temperatures. 
Middle: Rescaling of time to ct/ log t showing saturation at 
CT sat / log r sat ~ 1. Lines are the random walk prediction 
~ t/logt and a horizontal plateau showing saturation at 
C X2,2 ~ 1-8- Bottom: Rescaling of time to c 2 t/logt show- 
ing decrease of correlations which scales as the persistence 
function in Fig. 1111 



where D is the diffusion constant of the random walkers. 
The diffusion constant for dimers in the froth model will 
be of order unity since the moves carry no energy penalty. 
In Fig. E| we show reasonable fit with D = 1. We note 
that the dimers diffuse in a zig-zag fashion. A given 
dimer makes random steps of length unity along two non- 
orthogonal directions that depend on the relative position 
of the two defects that form the dimer (this orientation is 
constant as the dimer diffuses) . This zig-zag motion will 
lead to a reduction in the diffusion constant. Since there 
are only three possible dimer orientations, there is also 
the possibility of anisotropy in the correlations, but this 
was not observed in our simulations. Since the scaling of 
Fig. EDi s good and the fit reasonable, we do not attempt 
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FIG. 13: Scaling plot of (tC^ a/xi.a) as a function of the 
reduced variable r 2 /t. The full line is the independent random 
walk prediction 13111 derived in Ref. |27|. The fit is valid for 
r > 1 so the deviations at small r are greatest at small t. 
The deviations at large r are of the order of the numerical 
uncertainty in our simulations, but might also be an effect of 
the zig-zag motion of the dimers. 



any more rigorous analysis of these effects. 

Finally we consider the behaviour of the dynamic 
lcngthscale at saturation. The number of sites visited by 
a walker increases as (t/\ogt). Since the absorbing sites 
are uniformly distributed the walker will be absorbed on 
a timescale satisfying r sat ~ c _1 logr sat . In that time the 
walker travels an average distance 



.1/2 



-1/2 



[logc 



O (log log c- 1 )] 1 / 2 (32) 



In Fig. E| (top) we present scaling of the saturated cor- 
relation function according to this law. We also show 
in Fig. E] (bottom) that when times become of the or- 
der of or larger than the persistence time, the change in 
the susceptibility x'2 2 W i s mainly due to a change in the 
strength of the correlations while the dynamic lengthscale 
£2 2 i s approximately constant, just as in the plaquette 
model. We attribute the weak decrease of the correlation 
length at the longest times to the fact that the persistence 
function retains rather little information for times much 
greater than the relaxation times r. As a consequence 
the free random walk calculations in Ref. [27j do not ap- 
ply at large times because they do not take saturation 
effects into account. 

To summarise the results of this subsection, we have 
shown that the froth model behaves similarly to the pla- 
quette model with a mean free path for dimers that scales 
as cr 1 / 2 and a relaxation timescale that scales as c -2 with 
important logarithmic corrections due to the dimension- 
ality of the random walks performed by freely diffusing 
dimers. 



IV. CONCLUSION 

The plaquette model is a spin model with simple dy- 
namics and no finite temperature thermodynamic singu- 
larities. Its Hamiltonian gives rise to effective kinetic 
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FIG. 14: Top: Plot of C' 2 ^ m the state with a satu- 
rated correlation area showing good scaling as a function of 
r^/clog(l/c). Bottom: Plot of C^M a t logarithmically 
spaced times between 0.3r and lOr, with (5 = 7. As time 
increases, the strength of the correlations decreases by two or- 
ders of magnitude. On the other hand, the gradient of these 
traces at large distances gives the correlation length which 
is decreasing rather slowly with time (we show dashed lines 
marking e~ r ^ 25 and e~ r ^ 20 ). We constrast this behaviour with 
the unsaturated regime in which the correlations are always 
strong, and the correlation length increases as i 1//2 . 



constraints, and therefore to dynamical heterogeneity. 
Despite simple thermodynamics in the plaquette repre- 
sentation, spins have static correlations whose length- 
scales diverge at low temperature. The symmetries of 
the Hamiltonian mean that only rather specific corre- 
lation functions have non-zero expectations. However, 
measurement of the fluctuations in two-point quantities 
allow extraction of the relevant lengthscales. 

Spatial correlations are strongly anisotropic. For static 
correlations, this is the result of a symmetry of the spin 
system. Dynamical correlations are anisotropic because 
of a microscopic conservation law. We have identified 
three lengthscales in the plaquette model. The mean 



-1/2 



con- 



spacing between excited plaquettes, £| ' 
trols four point static correlations. Fluctuations in two- 
point structure factors are controlled by £;f at ~ c_1 
which represents the mean distance between adjacent ex- 
cited plaquettes in the same row or column. Finally there 
is a dynamical lengthscale £2,2(2) whose behaviour is dif- 
fusive at small times and saturates to £| tat at larger times. 
We also identified two timescales in the system. The 



first is the saturation time, r sa t 



which sepa- 
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rates the regime in which dynamical correlations have 
an increasing lcngthscale from the regime in which their 
lengthscale is saturated. The second is the relaxation 
time r ~ c~ 3 which represents the typical time for the 
spin field to lose the memory of its initial configuration. 
This separation of timescales arises from the non-trivial 
conserved quantities at the zero temperature dynamical 
fixed point. We also showed that these effects lead to 
unusual aging behaviour. 

We explained the dynamical length and timescales 
in the plaquette model in terms of diffusion of excita- 
tion pairs along one-dimensional paths. We showed that 
the froth model correlations are obtained by generalis- 
ing these results to allow the excitation pairs to diffuse 
in two dimensions. This introduces logarithmic correla- 
tions, but once these have been taken into account the 
scaling is as expected for non-interacting random walk- 
ers. 



The results connecting static and dynamic correlations 
in the plaquette model demonstrate one specific mecha- 
nism by which a 'mobility field' emerges from a more 
familiar interacting spin field. They are relevant to the 
description of structural glasses in these terms. The fact 
that all of the dynamical correlations can be explained 
in terms of independent random walkers provides further 
evidence that this behaviour is rather universal in sys- 
tems with dilute diffusing defects. 
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